Coulomb zero bias anomaly for fractal geometry and conductivity of granular systems 

near the percolation threshold. 
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A granular system slightly below the percolation threshold is a collection of finite metallic clusters, 
characterized by wide spectrum of sizes, resistances, and charging energies. Electrons hop from 
cluster to clusters via short insulating "links" of high resistance. At low temperatures all clusters 
are Coulomb blockaded and the dc-conductivity a is exponentially suppressed. At lowest T the 
leading transport mechanism is variable range cotunneling via largest (critical) clusters, leading 
to the modified Efros-Shklovsky law. At intermediate temperatures the principal suppression of a 
originates from the Coulomb zero bias anomaly occurring, when electron tunnels between adjacent 
large clusters with large resistances. Such clusters are essentially extended objects and their internal 
dynamics should be taken into account. In this regime the T-dependence of a is stretched exponential 
with a nontrivial index, expressed through the indices of percolation theory. Due to the fractal 
structure of large clusters the anomaly is strongly enhanced: it arises not only in low dimensions, 
but also in d = 3 case. 

PACS numbers: 72.23.Hk, 73.22.-f, 72.80.Tm 



Granular materials play important role in modern 
technology and material science (see [IJ). In the recent 
years nanocomposite granular materials were invented, 
with characteristic grain size a on the scale of 1 — 10 nm 
2]. For systems with so small grains quantum effects 
should be essential, in particular the Coulomb blockade 
effect Therefore it is important to understand, how 
the Coulomb blockade is manifested in realistic disor- 
dered granular metal. 

There are two large families of granular metallic mate- 
rials. In the systems of the first family conducting grains 
are randomly embedded in an insulating matrix (Fig[T^,), 
while the systems of the second family are mixtures of 
conducting and insulating grains (Figd})). The perco- 
lation 0, [i| is a general geometric phenomenon, generic 
for systems of both families. Some conducting grains may 
touch each other Q , establishing a good contact (with di- 
mensionless conductance G), while the conductances gij 
between grains i and j which do not touch each other, 
are much smaller: <C G. If there is a percolation 
via a network of touching each other conducting grains, 
then electrons can travel throughout the system hopping 
from grain to grain only via good contacts G. Otherwise 
hopping through some bad contacts g is unavoidable. 

If all conductances are small (both G <C 1 and g^ <c 
1), then the Coulomb blockade effect exists at each metal- 
lic grain of the system. It is characterized by charging 
energy on the scale Eq ~ e 2 /a. The mechanism of the 

electronic transport at T <C Eq in this case is either di- 
rect intergrain hopping, described by the Arrhcnius law 

a oc exp{~E act /T}, T ES « T « T Arr , (1) 

or variable range cotunneling (VRC)0, [H, 0], described 



by the modified Efros Shklovskii law [Toj ] 

aocexpj-^Es/T) 172 }, T«T ES , 



where 



rp t 1 s^ ) / r* 



C* ~ ln(l/7/), C(T) In [T| S /(T 2 + T ( 



(2) 

(3) 
(4) 

(5) 



Here g is a properly averaged intergrain conductance. 
The temperature T ei ~ [E^S] 1 / 2 < T ES [S being 
the typical level spacing in a grain) corresponds to the 
crossover from the inelastic cotunneling (at T > T e i) to 
the elastic cotunneling (at T < Tei). The details can be 
found in Thus, £ is a moderately large logarithmic 
factor and the Arrhenius law can only be observed in a 
restricted intermediate temperature range. In the low- 
G case the presence or absence of the percolation in the 
system is only relevant for the value of g, appearing in 
C* as an argument of the log-function, and, therefore, it 
is only of a secondary importance. 

In this paper we will study the large-G case: g^ <C 1, 
but G » 1. Here the mechanism of transport is very 
sensitive to the percolation transition. The resistance 
between two contacting metal grains is so low that the 
charge easily spreads over clusters of connected metal 
grains (see FigQJ, and it is clusters - not the individual 
grains - that may either be Coulomb blockaded, or not 
blockaded. If the fraction of the metal in the system 
x is larger than the percolation threshold x c , then the 
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FIG. 1: Clusters of conducting grains in granular materi- 
als, (a) - Metal grains (shown black) in an insulating matrix 
(shown grey). Randomly distributed distances dij between 
two non-touching conducting grains determine corresponding 
tunnel conductances gij. (b) - A mixture of conducting and 
insulating grains. The "links" (one-grain insulating bridges) 
between two conducting clusters are supposed to have the 
same conductances = g, while all longer bridges have 
gij -C g and can be neglected. 



Coulomb blockade of finite clusters is not relevant, since 
the infinite cluster of conducting grains exists (see [H, 
5]), the current goes through this infinite cluster, and 
electrons do not have to visit finite clusters whatsoever. 
It is not the case below the percolation threshold, for x < 
x c , where the infinite cluster does not exist, and electrons 
have to hop from one conducting cluster to another due 
to tunneling through high-resistance insulating bridges 
between them. 

Close to the percolation threshold the distribution 
N(n) of numbers of grains n in a cluster has a long power- 
law tail: N(n) ~ n~ r at 1 -C n <C n cr . This tail is cut 
off only at n ~ rt cr , where 



£, d f ~ (x c - x\ 



> 1 



(0) 



is the number of grains in a critical cluster, £ ~ (x c — x)~ u 
is its radius (measured in the units of a), df is the fractal 
dimension of the infinite percolation cluster (and of any 
large finite cluster with n in the range I -C n < n cr as 
well) . The values of relevant critical exponents are given 
in the table Q] An estimate for typical chargin g en ergy 
for a cluster consisting of n grains was found in [ll| : 



E c (n) 



E (o) n lv+.)/*d, 



We will see in what follows that the critical clusters with 
low charging energy 



E, 



(cr) 



E$\x c -x)^ 



(0) 

c 



(8) 



play the key role in low temperature transport at x < x c . 

It is natural to expect that the universal character of 
the percolation transition, generic for most granular sys- 
tems, should lead to the essential universality of the con- 
ductivity mechanism for i c - i < 1. One should, how- 
ever, have in mind that besides the topological disorder 
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TABLE I: Numerical values of some critical exponents 



(that is responsible for the percolation phenomena) , there 
is yet another disorder: the randomness of "bad" conduc- 
tances g^. Indeed, oc exp(— 2Kdij) where dij is the 
separation between the two grains and k is the tunneling 
decrement of the electronic wave function in the insula- 
tor. If typical value of redy 3> 1, then the dispersion of g 
is exponentially wide. 

In mixtures the distribution of g is approximately dis- 
crete: the thickness d.^ of the insulating interval between 
two metal clusters is, roughly, measured in the units of 
the diameter of the insulating grain d. Then the tunnel- 
ing conductances of shortest one-grain insulating bridges 
(links) are g\ oc exp(— 2nd)] the conductances of two- 
grain bridges are gi oc exp(— And) <C gi, and so on. For 
such a system one can simply ignore all long insulating 
bridges and take into account only the shortest - the 
links, ascribing the same conductance g = g± to all of 
them. The conduction process in the resulting system 
resembles the next nearest neighbor (NNN) percolation 
(see ll| for detailed discussion of this process). 

For metal grains, embedded in an insulating con- 
tinuum, the distribution of dij is essentially continu- 
ous. This fact introduces to the system yet another 
percolation-like physics, similar to that of the standard 
hopping conductivity (see, e.g., [12]). In the present pa- 
per we do not consider this facet of the problem explic- 
itly, so that, strictly speaking, the consideration below 
is directly applicable only to mixtures. We expect, how- 
ever, that the principal universal features of conductivity, 
based on the universal properties of the clusters distri- 
bution near the threshold, will be present also for the 
systems with continuous insulating matrix. 

The temperature dependence of conductivity of a mix- 
ture in the intermediate range of temperatures 



E. 



(cr) 



e { ™ ] «T« J eg )) , 



(9) 



(7) was already studied in the previous paper |ll| . Here 



is the charging energy of a "marginal cluster" , whose 
classic resistance R m ~ 1 . The classic resistance across a 
fractal cluster of n grains (see FigfJJand Refs. [1, EH) is 

R(n) - G-V^ 2 -^/ 1 ^, (11) 

so that the number of grains in the marginal cluster is 

n m - G" d f /[fl+i2 - d),/] > 1. (12) 
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FIG. 2: Definition of the classic resistance R across a cluster. 
External voltage V is applied to points A and B at the pe- 
riphery of the cluster, then the ratio of V and the current I 
in the external circuit is the resistance R. 



The clusters with n < n m (and, therefore, with R < 1) 
can be treated as structureless point-like supergrains 
characterized by the unique quantity - the charging en- 
ergy Ec(n). Under the condition (J9j) that is true for all 
relevant clusters in the system and the conductivity is 
described (see [ll]) by the formulas 



c 
c 0) 



T < T croS8 (a;), (13) 
), (14) 



where 



4 0) K*c 



■ x)/A c 



i(a»+s)/o 



(15) 



is the temperature of the crossover from the "insulator- 
controlled" conduction mode, where the resistivity is 
dominated by the insulating links, to the "critical 
crossover" mode, where the voltage drops occur both on 
links and on the conducting clusters. Note, that in the 
latter regime the conductivity cr C ross does not depend on 
x. The critical exponents 



Q = (d-2)v + s-l Q ,_ /ie 



(16) 



(0) 



are given in the table Q] At high temperatures T ~ E\ 
when the Coulomb blockade effect becomes irrelevant, 
the expressions (|13|14[) match with the known results (see 
[HI]) obtained in the absence of the Coulomb effects: 



a; 



ins 
(0) 



g{x c 



(A cross «i c -i«1), (17) 



( Xc - x <A cmss ), (18) 



where the width of the critical crossover domain of con- 
centrations is 



(g/G) 



(19) 



The physics behind the results (|13I14[) is as follows: In 
the temperature range ([9]) it is possible to find conduct- 
ing paths consisting of only large metal clusters connected 
by links, so that an electron never visits the Coulomb- 
blockaded small clusters with n < ucb(T) (and with 
charging energies Ec{n) > T). Provided nc^,(T) < n a , 
the NNN-percolation in the system persists despite the 
fact that all small clusters with n < uqb (T) are not avail- 
able for travelling electrons. This fact is by no means 
trivial, since almost all conducting grains in the system 
belong to small clusters. The power-law suppression fac- 
tor (T/E^) e in the conductivity reflects just the reduc- 
tion of the number of available conducting paths. 

In the present paper we will address the following ques- 
tions: 

1. What is changed if T <C Eq, so that the rel- 
evant clusters are not point-like zero-dimensional, 
but essentially extended objects with their internal 
dynamics? How one should describe the Coulomb 
blockade phenomena in this case? 

2. What is the conduction mechanism and what is the 
T-dependence of conductivity at very low temper- 
atures, when all clusters are Coulomb blockaded? 

The low-temperature tunneling into a finite, but ex- 
tended conductor is suppressed due to the process of 
charge spreading which transforms the initial point-like 
distribution of the tunneling charge into the smooth 
equipotential distribution. The corresponding sup- 
pression factor exp(— S spr (T)) is nothing else, but the 
Coulomb zero-bias-anomaly (ZBA) factor, which appears 
in the tunneling probability alongside with the usual 
Coulomb blockade factor exp(— Ec/T). At low tempera- 
tures, when the sample size L is smaller than the "spread- 
ing length" L(T), the activation factor dominates, while 
at intermediate temperatures (for L > L(T)) the prin- 
cipal contribution comes from the ZBA factor. Explicit 
expressions for S spT (T) (as well as for L{T)) are known for 
diffusive conductors of different geometries: one-, two-, 
and three-dimensional (see [3 H H E3, H 0, ) ; 
the suppression is exponentially strong only in Id case. 
In our problem we deal with an unusual case: we have to 
find the ZBA-factor for tunneling between two large clus- 
ters with fractal geometry, where the diffusion constant 
D, the conductivity a, and the dielectric constant e are 
scale-dependent. We will see that for such a geometry the 
suppression is exponentially strong in any dimension: in 
particular, in 2d and in 3d. 

Let us consider a tunneling of an electron through a 
link between two clusters with similar numbers of grains 
n% ~ ri2 ~ n. Immediately after the intercluster hop 
the system finds itself under the barrier, with large en- 
ergy deficit ~ Eq^ , so that the charge density has yet to 
spread over a large region in a tunneling manner, before 
the system manages to get from under the barrier (see 



4 




FIG. 3: The charge spreading process. After the tunneling 
of an electron between two clusters (black and grey) through 
one of the links (marked by a cross) a diffusional spreading of 
charge occurs in both clusters. 



FigGU). As a result, the suppression factor has the form 
exp{-S(T)} with 



S(T,n) = E c (n)/T + S spr (T, n), 



(20) 



where the under-barrier action S spr (T,n) can be esti- 
mated by the semiclassical method. Different variants 
of such method were proposed in , Il7| , in this pap er 
we will use the variant due to Levitov and Shytov [17| . 

Strictly speaking, the act of the electron hop between 
two neighboring large clusters creates an electron-hole 
pair, and each component of this pair then spreads 
over its own cluster. However, this separation does not 
prevent the Coulomb interaction of both components. 
The electron-hole interaction effectively leads to partial 
screening of self-interaction in each cloud. If both clouds 
develop in the same spatial domain, then this screening 
is strong and leads to a dramatic suppression of the ac- 
tion S spT . Exactly this situation arises, when electron 
and hole clouds proliferate in nearby parallel planes with 
similar properties. The action in this case is parametri- 
cally smaller than that in the STM case, where the hole 
is immediately evacuated [17 1. 

In our problem we apparently have an intermediate 
case. The neighboring large clusters to some extent in- 
terpenetrate, but by no means they coincide. The rela- 
tive overlap is probably considerable, but definitely it is 
not close to a complete coincidence. Obviously, the effect 
of mutual screening is important for our problem, but it 
only leads to an effective suppression of self-interaction 
by a numerical factor, without introducing any new scale 
or any new small parameter. Thus, we conclude, that 
neglecting the mutual screening would overestimate the 
under-barrier action S spl only by a numerical factor of 
order, say, two. Since we are anyway not able to deter- 
mine the numerical factor in the action, we will neglect 
the effect of mutual screening in what follows. 

We write the expression of the action S spl in the case, 
when the process is controlled by finite temperature (not 
by the external voltage), and neglect the mutual electron- 



hole screening: 

S spl (T,n) 
d d q 



2ttT 



^ Q 2irT(2k+l) + D q q 2 



i /L (2Tr) d 2TrT(2k+l) + a q q 2 U q 



D q q 2 



(21) 



(see [20j for details). Here L is the size of an electrode, 
a and D are its conductivity and diffusion constant, and 
U is the screened Coulomb interaction. The wave-vector 
q is measured in the units of 1/a. Now we extend the 
result (|2T1) to the fractal case with 



Gq 

U a 



D q ~ GE^q* 



E (0) 1-d+s/u 



L 



(22) 



(0) 



4 W)- 2 <<4 0) , 

therefore the diffusion terms D q q 2 can be neglected in 



(see [j|). The energy scale E 



both denominators in (|2"Tj) . The sum over k in (|2~4"|) is 
dominated by k ~ 1 and the integrals over q are domi- 
nated by q ~ i(T) _1 , where 



L(T) ~ (gE^/t) 



v/[n+s+{Z-d)v] 



(23) 



The spreading factor is the dominant one, if 1 <SC 
L(T) < L. Introducing new variable Q = qL(T)(2k + 
l^-v/[v+s+(3-d)v] ^ anc j extending the integration to the 

infinity, we get 



S spr (T,n)~(4 m) /T)' 



(24) 



where the new critical exponent tp is expressed in terms 
of the universal indices of the percolation theory: 



fX+(2- d)v 
H + s + (3 — d)v 



0.33, 
0.41, 



(d = 2), 
(d = 3). 



(25) 



The situation, when the action (|20|) is dominated by 
the activation term, we call the "Coulomb blockade" , 
while the regime, dominated by the spreading term is 
called the "Coulomb ZBA". The phase diagram for dif- 
ferent regimes on the n — T plane is shown in FigHl The 
crossover line is defined by 



ncB (T) 



(T/E ( c ] )- udt/{v+s \ (T > ^ m) ), 
(T/GE^y^f^+s+V-W, (T « E^). 



The crossover temperature Tai-i-Mi below that all clus- 
ters in the system are Coulomb blockaded, so that the 
Arrhcnius activated regime of conduction sets on, can be 
found from the condition n cr (x) = ncB(7Arr): 



T A rr{x) 



E. 



(0) 



G(x c 



{x c 



x > A m ) 



x < A m ), 



(26) 
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No Clusters 




FIG. 4: Different regimes for tunneling between two clusters 
with n grains each on the n — T plane. The most interesting 
case n cr > rim is shown; in the opposite case, for n cr < n m , 
there are no clusters in the "Coulomb ZBA" domain. 



where 



A m = 



G -l/[^+(2-dH 



(27) 



and x m is defined by the equation n cr (x m ) = n m . 

The result (f2~4")) is very different from the standard 
Coulomb ZBA effect in plain (non- fractal) systems, 
where the stretched exponential law (with power 1/2 ) 
appears only in one-dimensional case B m M m, 
while in 2-d there is only a log-squared function of T in 
the exponent and in 3-d the anomaly is weak. The rea- 
son for such a striking difference is the nontrivial scale 
dependence (|22| of both dielectric screening and the con- 
ductivity in the fractal system. 

In the range T Arl (x) < T < E^ n) the NNN- 
percolation via large non-blockaded clusters (with 
n CB(T) < n < n cr ), is still possible. However, the tun- 
neling between these clusters is strongly modified due to 
the Coulomb ZBA. As a consequence, the expressions for 
conductivity (|13I14|) are modified also: The result ^ 
for the insulator controlled transport is changed to 



CTi, 



(x,T) 



T 



E, 



(0) 

c . 



B exp 



E. 



(m)' 



T 



(28) 



where c is some unknown numerical constant, B = 
B \ T/Eq \ is unknown prcfactor. Note that the ex- 
pression on the right hand side of (f2"8"|) does not depend 
on x. The result (fT4")) for the conductivity in the critical 
crossover range is changed to 



T 



e' 



changed to 



a (0) 

u cross 



E, 



(o) 



exp 




and the expression for the crossover temperature (| 1 5[) is 



3 (:r)~4 m) ln 1 /*' [A* CIO J(x c - x)} 



The modified critical exponents 



■ ( fl + s)lM + (2-d)^] 



e = 



(d-2)v + s-l 
fi + s + (3 — d)i/' 



e' = 



/' 



(30) 
(31) 



(32) 



Now we turn to the question about the nature and 
the temperature dependence of the conductivity at low 
temperatures T < T^r- Here an electron can not avoid 
Coulomb blockaded clusters, since practically all clus- 
ters are blockaded. At moderately low T electrons travel 
through the "critical network of critical clusters" . This 
mode of transport is similar to the standard nearest- 
neighbor-hopping regime in the hopping conductivity 
(see [Hj])- The hopping goes via critical clusters, since, 
on one hand, they have small charging energies ~ E c 



(cr) 



and, on the other hand, they form an NNN-percolation 
network (i.e., it is possible to travel by hopping from one 
critical cluster to another through direct links). In this 
regime the conductivity obeys the Arrhenius law |T]) with 



E, 



E, 



(cr) 



act {%) ~ 

If the temperature is lowered furtherjthen, in the spirit 
of the Mott's variable range hopping |12L|2l| the variable 
range cotunneling regime (VRC) @, H, Q sets on. This 
regime is characterized by distant hops between resonant 
clusters with small charging energies. The hops are real- 
ized as acts of multiple cotunneling (elastic or inelastic) 
through chains of nonresonant clusters, connecting the 
resonant ones. Thus, VRC involves two sorts of clusters: 

• The terminal ones, where the real charged states 
can occur. These clusters should be resonant. 

• The intermediate clusters, forming chains of vir- 
tual intermediate states. These nonresonant clus- 
ters should be as large as possible - to minimize the 
number of intermediate states. Of course, they also 
should be NNN-connected. 

What is the nature of the resonant clusters? While the 
characteristic scale of the charging energy for a cluster 
of size n is Ec(n), its particular value for a given clus- 
ter is random, due to random form of the clusters, and, 
most important, due to random electrostatic potentials of 
the surrounding. The electrostatic interaction of metallic 
clusters with charges Qi (measured in the units of elec- 
tronic charge e) is H C {Q} = %■ J2ij [ C ^] ij QiQj, where 
Qi = (Qi — qi) and CV, is the matrix of capacitances. The 
so called "background charges" —1/2 < qi < 1/2 are due 
to random potentials of stray charges, trapped in the in- 
sulator; the set of qi is chosen so that the ground state 



6 



corresponds to Qi = 0. When one extra electron (one ex- 
tra hole) is put on the cluster i, the energy of the system 
is changed by 



E, 



(±) _ 
Ci ~ 



(33) 



The resonant clusters are characterized by anomalously 
small Eq) or E^. Counter-intuitively, the principal 
contribution to the low energy density of states comes 
not from large resonant clusters, but from small reso- 
nant clusters. To prove this, let us first neglect the long 
range Coulomb interaction (i.e., the Coulomb gap effect) 
for a while, and consider the so called "density of ground 
states" Vqs{E) (see Q for the detailed definition). It 
is instructive to decompose it into the sum of contri- 
butions VQ S (E,n) of clusters with fixed n. Under the 
most natural condition of strong charge disorder (when 
qi are homogeneously distributed in the interval — 1/2 < 
qi < 1/2) each i/qs(E, n) is a structureless function with 
a single scale Ec(ri), so that v^ s {E,n) ~ 1/Ec{n) for 
E < E c (n) and v { °l(E,n) « for E > E c {n). As a 
result, for the sum we obtain 



n(E) 



« J2 n- T /E c {n) « ^ n -[r-(s+u)/ud f ] _ L ^ 

n— 1 n 

Here ri(A) is defined by Ec{n) = E. Since r — (s + 
v)/vdf > 1, the sum is indeed dominated by small clus- 
ters with n ~ 1. It is important to stress that, although 
most clusters with low charging energies are small reso- 
nant clusters, the latter do not play any role in the con- 
duction processes at T > Tj^ IT . Unlike large clusters, 
these small resonant clusters are not NNN-connected: 
they are separated from each other by many nonresonant 
ones. That is why we did not take them into account in 
the previous parts of this paper (as well as in 

Taking into account the long range Coulomb interac- 
tion modifies the result (|34p and leads to the appearance 
of the soft Coulomb gap [71 Il2|. 



v GS {E) ~ \E\ d -\e/e 



(35) 



but it can not alter the small-cluster nature of the ma- 
jority of low energy states. 

The asymptotic law ([33]) is valid only for lowest ener- 



gies \E\ < E { "\ while in the range Eg* 1 < E < 1 the 
density of states is strongly modified by the fractal char- 
acter of the system. However, the VRC regime is actual 
just in the temperature range, where the resonant clus- 
ters have energies E -C Eq r ' , so that the formula is 
sufficient for our purposes. 



(cr) 



The intermediate clusters involved in the acts of mul- 
tiple cotunneling are the critical ones. In contrast with 
small resonant clusters, which are always point-like, these 
clusters can be either point-like (if n CI < n m ), or ex- 
tended (if n m < n cr ). There is a considerable difference 
in the cotunneling probability between the two cases. 

Under the condition n cr < n m all the critical clusters 
at low temperature act as effective " supergrains" with 
characteristic charging energy E^\ The average num- 
ber of links, connecting two adjacent critical clusters, is 
Minks(C) ~ VO^c — x) (see [n|). Therefore the effective 
conductance between the two clusters is 



9cS ~ Minks (£) ~ 9 1 (x c - x). 



(36) 



In this paper we consider only the case g e s <C 1. Wc 
also concentrate on the case the inelastic cotunneling; a 
subtle question about the elastic cotunneling and the x- 
dependence of T e i(x) will be discussed in a separate pub- 
lication. For point-like clusters the method of evaluation 
of the probability of multiple cotunneling, proposed in 
[8|, can be applied directly. Introducing the characteris- 
tic number K of critical clusters, separating two resonant 
ones, and the width A of the Mott strip, we obtain 



fins oc g% 



(A/A' 



E, 



(cr) 



2 A' 



exp{-A/T}, 



(37) 



for T e i(x) <C T <C TArr(^)- Then, having in mind the 
relation A • ^Gs(A)(^A) d ~ 1 between A and A, and 
optimizing the conductivity (|3T|) . we arrive at the formula 
© with A ES ~ C{T)E { ™ ] and 

C(T) =C*+2 In {T ES /T) , C* - In (l/ ffoff ) . (38) 

Thus, for point-like critical clusters Tes ~ C* . 

In the case n cr < n m we have to deal with cotunnel- 
ing through extended objects with large resistance. Such 
cotunneling was considered in [2(| for the case of long 
diffusive wire; the corresponding exponential factor was 
found by means of modified Levitov-Shytov scmiclassi- 
cal method. The approach of [2(| can be applied also to 
the extended fractal clusters. Writing the action for the 
process, where an electron and a hole are simultaneously 
injected into an extended cluster at its opposite ends, at 
distance r ~ £ from each other, we obtain 



Sr. 



2ttT 



d d q 



" Q 27rT(2fc + 1) + D q q 2 
L/ g sin 2 (q.r/2) 



i/ 4 (2^) d 2vrT(2fc + 1) + a q q 2 U q + D q q 2 

~ 2R(n cr ) In (T Arr /T) . (39) 

The expression (|39|) differs from ((2T|) only in the factor 
sin 2 (q-r/2). In the "Coulomb ZBA scenario" the integral 
over q in (|2"Tj) is cut off at the intrinsic scale q ~ A _1 (T), 
controlled by the temperature and a:- independent. On 
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the contrary, in the VRC regime L(T) £, so that the 
integral is cut off at q ~ £~ 1 . As long as we are interested 
only in the inelastic cotunneling, the diffusion terms D q q 2 
again can be neglected in both denominators in (|39|) 

Now, repeating the arguments of the formula ([57]) 
with account for the exponential suppression of the co- 
tunneling amplitude, we get again the result (J2J), but 
with E ES ~ C{T)E ( " ] R{n„), Tes ~ T A „/£* and C* ~ 
l/i?(n cr )ln(l/<?ofr). 

Finally the results for both n cr < n m and n cr < n m can 
be unified in the following form: the crossover tempera- 
tures Ta it and Tes, separaiting the domains of validity of 
the Arrenius law |T]) and the Efros-Shklovskii law ([2]) are 
strongly reduced in the vicinity of percolation threshold, 



CLASSIC (1) 



CLASSIC (2) 



Tat 



E. 



C f-tUiS, 



(40) 



C* = 1 + (l/i? eff ) In (l/ ffe ff) , i?eff « 1 + i?(n cr ), (41) 



#(n cr ) - G _1 (a; c - x )-Vl»+(2-d)»] ^ 



(42) 



and i? c ff has the meaning of the effective dimensionless 
resistance across the critical cluster. 

Note that the gap between Ta T t an d Tes only exists, 
if £.* » 1. This condition is fulfilled for x c — x > A , 
where 



A = A m [Hl/g)]-^ 2 - d ^ 



(43) 



while for x c — x < A C* ~ 1 so that the gap shrinks 
to zero and there is no room for the Arrhenius law (see 
Fi^SJ). The "activation energies" E act and Ees entering 
([!]) and are also reduced: 



E act ~ E^\ E ES ~ £(T)E ( ">R eS , (44) 



(«); 



£(T) = £* + 2 In (T ES /T), 



(45) 



In conclusion, we have studied the low-temperature be- 
havior of a granular system near the percolation thresh- 
old. The peculiarities of this behavior stem from the 
strong dispersion and fractal properties of the conduct- 
ing clusters. As in the system away from the percolation 
threshold, the transport mechanism at lowest tempera- 
tures is variable range cotunneling, and the correspond- 
ing temperature dependence of a is the Efros-Shklovsky 
law. However, the parameters entering this law are dra- 
matically renormalized in the vicinity of the threshold. In 
particular, the onset of the VRC regime is shifted to lower 
temperatures. At higher temperatures the VRC is re- 
placed by alternative phisical mechanisms. One of them 
- the NNN-percolation with excluded small (Coulomb 
blockaded) clusters, was studied in [llj, the other - the 
Coulomb zero-bias anomaly scenario is first discussed in 
the present paper. The Coulomb interaction suppresses 



\ "EXCLUDED 
"EXCLUDED CLUSTERS" (1) J* \ CLUSTERS" (2) 




Across (A cr0 ss) (x c -x) 



FIG. 5: Phase diagram of different regimes for the conduc- 
tivity. "Classic (1)" - classic T-independent transport, con- 
trolled by the links (see EqJTTJl ) ; "Classic (2)" - the same in 
the critical crossover region (see Eq{TS}); "Excluded clusters 
(1)" - transport via large (non-blockaded) NNN-connected 
clusters, the T-dependence is a power law (see Eq (|13j) '); "Ex- 
cluded clusters (2)" - the same in the critical crossover re- 
gion (see Eq (|14p h "Coulomb ZBA (1)" - transport is con- 
trolled by the Coulomb zero-bias anomaly, the T-dependence 
is stretched exponential (see Eq (|28[l ); "ZBA (2)" - the same 
in the critical crossover region (see Eq (|29p ); "Arrhenius" - 
direct hopping between Coulomb blockaded critical clusters, 
the T-dependence is simple activation (see Eq|T}); "VRC" - 
variable range cotunneling (elastic or inelastic) between small 
resonant clusters via chains of critical non-resonant clusters, 
the T-dependence is the modified Efros-Shklovskii law (see 
Eq©). 



the probability of tunneling between large metallic clus- 
ters in the "zero bias anomaly" manner. Due to the frac- 
tal structure of clusters the temperature dependence of 
the Coulomb ZBA factor is described by stretched expo- 
nential law (|28p with nontrivial index ip. 

The Author is indebted o M.V.Feigelman for numerous 
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cial for this work. This research was supported by RFBR 
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